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ABSTRACT 

We present measurements of the optical and X-ray continua of 108 AGN (Seyfert 
Is and quasars) from the Rosat International X-ray/Optical Survey (RIXOS). The 
sample covers a wide range in redshift (0 < z < 3.3), in X-ray spectral slope 
(-1.5<a x <2.6) and in optical-to- X-ray ratio (0.4<a ox <1.5). A correlation is found 
between a x and a ox ; similar correlations have recently been reported in other X-ray 
and optical samples. We also identify previously unreported relationships between the 
optical slope (a op t) and cv x (particularly at high rcdshifts) and between a op t and a ox . 
These trends show that the overall optical-to-X-ray continuum changes from convex 
to concave as a x hardens, demonstrating a strong behavioural link between the op- 
tical/UV big blue bump (BBB) and the soft X-ray excess, which is consistent with 
them being part of the same spectral component. 

By constructing models of the optical-to-X-ray continuum, we demonstrate that 
the observed correlations are consistent with an intrinsic spectrum which is absorbed 
through different amounts of cold gas and dust. The intrinsic spectrum is the sum of 
an optical-to-soft X-ray 'big bump' component and an a x —l power law; the column 
density of the cold gas ranges from to ~ 4 x 10 21 cm -2 and the dust-to-gas ratio 
is assumed to be Galactic. The 'big bump' may be represented by a Tb rem ~ 10 6 K 
thermal bremsstrahlung or an accretion disk with a surrounding hot corona. The 
scatter in the data can accommodate a wide range in big bump temperature (or black 
hole mass) and strength. A source for the absorbing gas may be the dusty, molecular 
torus which lies beyond the broad line-emitting regions, although with a much lower 
column density than observed in Seyfert 2 galaxies. Alternatively, it may be the bulge 
of a spiral host galaxy or an elliptical host galaxy. 

Key words: Quasars: general - Galaxies: Active - Galaxies: Seyfert - X-rays: gen- 
eral. 



1 INTRODUCTION 

Three major continuum components are observed in the 
optical-to-X-ray spectra of non-blazar AGN; the big blue 
bump (BBB), which rises through the optical and UV 
(e.g. Edelson & Malkan 1986; Elvis et al. 1994), the 'canon- 
ical' hard X-ray power-law, which dominates above ~2 keV 
(e.g. Mushotzky 1984; Turner & Pounds 1989; Comastri 
et al. 1992) and the soft X-ray (E <1.0 keV) 'excess' above 
the X-ray power-law [e.g. Arnaud et al. 1985; Turner & 
Pounds 1989; Walter & Fink 1993 (hereafter WF)]. The 
BBB generally dominates the bolometric luminosity of non- 
blazar AGN spectra and the soft X-ray excess may be its 
high energy tail. The variability of the BBB indicates that 
it is emitted from the nuclear regions of AGN (e.g. Clavel 
et al. 1991) and a popular model is that of blackbody-like 
emission from an optically-thick accretion disc (AD; e.g. Sun 
& Malkan 1989; Czerny & Elvis 1987; Madau 1988). 



Measurements of the BBB's spectral form, and conse- 
quently model development, are limited however because 
the peak of the BBB, and much of its frequency coverage 
are lost in the EUV where only a handful of AGN may be 
detected. Nonetheless, using UV and soft X-ray data from 
X-ray bright Seyferts, WF reported a correlation between 
the 1375 A-to-2 keV flux ratio and the slope of the soft X- 
ray (0.2-2.0 keV) spectrum, in the sense that AGN with rela- 
tively soft X-ray spectra also have a high UV-to-X-ray ratio. 
They suggested that the correlation was due to changes in 
the flux of a single optical-to-soft X-ray 'big bump' relative 
to an underlying continuum. Laor et al. (1994; hereafter 
L94) have reported a similar correlation in their sample of 
UV-excess quasars. This has provided the first direct evi- 
dence of a link between the BBB and the soft X-ray excess, 
although there still remains doubt as to whether they are the 
same component, or separate but related (Fiore et al. 1995). 

One restriction of the three samples mentioned previ- 
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ously (i.e. WF, L94 and Fiore et al. 1995) is that they are 
dominated by low-redshift (z < 0.4) AGN whose proper- 
ties may differ from relatively high-luminosity quasars. The 
L94 sample is made up of nearby (z < 0.4) quasars with 
strong, blue optical/UV continua and the Fiore AGN have 
a relatively high optical-to-hard- X-ray (2500 A-to-2 keV) lu- 
minosity ratio; these criteria favour AGN with strong BBBs. 
The WF sample contains 58 soft X-ray (0.08-2.4 keV) bright 
Seyferts which may also favour AGN with strong BBBs. A 
sample of AGN which spans a wider range in X-ray slope 
and optical to X-ray ratio would be better suited to a more 
thorough investigation of the soft X-ray excess and the op- 
tical/UV BBB. 

We have used AGN (i.e. the Seyfert Is and quasars) 
detected as part of the RIXOS survey (Mason et al. 1996) 
to extend this research. The RIXOS survey, which is de- 
scribed in Section 2, is made up of serendipitous sources de- 
tected in long-exposure pointed observations made with the 
ROSAT Position Sensitive Proportional Counter (PSPC; Pf- 
effermann et al. 1986). RIXOS sources are selected accord- 
ing to their flux above 0.4 keV, thereby reducing the prefer- 
ence for AGN with a strong soft X-ray excess. They cover 
a wide range in redshift (0 < z < 3.3), X-ray spectral slope 
(-1.5<a x <2.6) and in the optical-to-X-ray flux ratio 
(0.4<a ox <1.5), providing an extensive coverage of param- 
eter space in which to search for trends and correlations. 

In Section 2 we outline the essential details of the sur- 
vey and the process for extracting and measuring the X-ray 
and optical data. The overall properties of the sample, cor- 
relations and 'mean' optical-to-X-ray spectra are presented 
in Section 3 while in Section 4 we present a model for the 
RIXOS AGN continua. The implications for the big bump 
and for the structure and geometry of AGN are discussed in 
Section 5. 

In total, 108 AGN are included in this analysis; the full 
source list giving X-ray and optical positions, fluxes etc., 
will be presented in a future paper (Mason et al. 1996) 
when the survey is complete. A companion paper to this 
(Puchnarewicz et al. 1996) will discuss the relationships 
between the optical and UV lines and the continuum pa- 
rameters. 

2 THE RIXOS SURVEY 

The RIXOS survey is compiled from serendipitous sources 
detected in pointed PSPC observations. To be included in 
the survey, an observation must have an exposure time of at 
least 8 ksec and it must have been taken at a Galactic lati- 
tude greater than 28° (where the Galactic absorbing column 
density is low). Only sources within 17 arcmin of the centre 
of the field and with a flux greater than 3xl0~ 14 erg sec -1 
in the ROSAT 'hard' band (0.4-2.0 keV) are used. These 
criteria ensure that source selection is made irrespective of 
the strength of any soft X-ray excess which dominates at 
low energies (E < 0.5 keV), however by ignoring any counts 
below 0.4 keV, there will be a bias against ultrasoft X-ray 
AGN like E 1346+266 and RE J1034+396 (Puchnarewicz, 
Mason & Cordova 1994; Puchnarewicz et al. 1995), which 
have relatively little flux in the 'hard' band. Most sources 
lie behind relatively low Galactic Nh (A^HGai), 60 per cent 
of the RIXOS AGN have an Af H Gai< 2 x 10 20 cm" 2 and all 
have an A?HGai below 5 x 10 20 cm -2 . 
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Figure 1. The distribution of total source counts for the spectra 
used in this analysis 

The RIXOS survey is 100% complete over 20 sq degrees 
to a flux of 8.4xl0~ 14 erg cm~ 2 s _1 and 93% complete over 
15 sq degrees to a flux of 3xl0~ 14 erg cm' 2 s _1 (Mason 
et al. 1996). The subsample analysed for this paper is drawn 
from the first objects to be identified from the whole sur- 
vey: confused sources (i.e. those which were less than 30 
arcsec from other RIXOS sources), sources with poorly con- 
strained fits to the X-ray data and objects with very weak 
optical spectra have been discarded. We find no significant 
differences between the X-ray parameters of the AGN pre- 
sented here and those from the 100% complete sample (see 
Mittaz et al. 1996), thus we expect that the properties of 
this subsample are representative of the final survey. 

2.1 X-ray data 

To measure the X-ray spectral distribution of each source, 
the PSPC counts are divided into three bands, S (0.1 to 
0.4 keV; channels 8 to 41), HI (0.5 to 0.9 keV; channels 52 
to 90) and H2 (0.9 to 2.0 keV; channels 91 to 201). At the 
first stage in the reduction process, periods of high count 
rate in the time series of the background were removed. The 
source signal in each band was measured by summing the 
counts in a circle centred on the source's centroid position: 
in most cases the radius of the extraction circle was 54 arc- 
sec, however where other sources lay closer than 54 arcsec, 
a circle with a radius of one-half of the distance between 
the source and its nearest neighbour was used (in these 
cases, extraction radii are indicated in column 2 of Table 
1). The 54 arcsec radius was chosen by simulation to opti- 
mize the signal-to-background ratio for the generally faint 
sources in the RIXOS survey. All source fluxes were cor- 
rected for counts falling outside the extraction circle using 
the formulae in Hasinger et al. (1994), out to a radius of 
3 arcmin. A large, annular region, centred on the point- 
ing position of the field, was used to calculate a model for 
the background, after correcting for vignetting and having 
masked out any source emission. 
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The broadband measurements were combined to pro- 
duce 'spectra' with three data points for each source. A 
histogram illustrating the distribution of the total counts 
in the sources used in the present analysis is plotted in Fig 

1 and shows that one half have at least 150 counts. The 
'3-colour' X-ray spectra were fitted using single power-law 
models, fixing the absorption to the Galactic column den- 
sity [this was calculated by interpolating between the 21 cm 
measurements of Stark et al. (1992)]. With two free pa- 
rameters (power- law index and normalization), this left one 
degree of freedom for the fits. 

Many of the RIXOS sources have less than 15 counts 
in at least one of the S, HI or H2 bands; at this level, 
the assumption of a Gaussian probability distribution for 
the data is inappropriate and Poissonian statistics must be 
used. Thus instead of finding the best-fitting model for each 
source by minimizing \ 2 (x 2 ls normally used when fitting 
X-ray data with a higher signal-to-noise ratio and assumes 
Gaussian statistics), the best-fit parameters (normalization 
and index) were calculated by minimizing the Cash statis- 
tic (Cash 1979) which is appropriate when considering data 
with relatively few counts. The errors were calculated from 
the 90 percent confidence contours on the fit, using a method 
identical to that for \ 2 (see Cash 1979). A model for the 
total number of counts in each band, which included the 
model of the background, was constructed for each source 
and compared with the actual data (i.e. source plus back- 
ground). The detector response function and corrections 
for vignetting effects and particle contamination were folded 
into the model. The results of this spectral fitting proce- 
dure have been verified by extensive simulation. A com- 
plete discussion of the X-ray spectral fitting procedure for 
AGN and other sources in the RIXOS survey (i.e. including 
AGN, emission line galaxies, clusters and stars) is presented 
in Mittaz et al. (1996). 

The best-fitting X-ray spectral indices and their 90% 
(Cash) errors are given in Table 1; a x , and indeed all spectral 
indices (a) are defined throughout such that F„ oc u~ a . 

2.1.1 Luminosities 

The logarithms of the luminosities at 0.2 keV (Lo.2kev) an d 

2 keV (I/2kcv) in erg s -1 Hz -1 in the rest-frame of the quasar 
were calculated using the best-fit power-law model for each 
individual source, assuming a value of 50 km s _1 Mpc -1 for 
the Hubble constant (Ho) and for the deceleration param- 
eter (go; these values are assumed throughout this paper). 
Errors on the X-ray luminosities have been calculated from 
the Cash 90% limits on the best-fitting power-law models 
(including uncertainties on both the slope and normaliza- 
tion); i .2kev, i2kcv and their respective errors are listed in 
Table 1. 



2.1.2 Spurious correlations 

Because the RIXOS sample is flux-limited, spurious correla- 
tions may be introduced between the various continuum pa- 
rameters (slopes and luminosities) and these must be borne 
in mind when considering the viability of trends and rela- 
tionships in the data. Sources are selected according to their 
flux in the hard (0.4-2.0 keV) band in the observer frame, 
thus 2 keV in the rest-frame of the AGN is contained within 
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Figure 2. The X-ray luminosities at (a) 2.0 keV (-L2kcv) and (b) 
0.2 keV (io.2kev) plotted as a function of X-ray spectral index, 
a x (filled circles). Also plotted (as open squares) are compari- 
son, random datasets (see Section 2.1.2) to illustrate the spurious 
correlation between Lo.2keV an d a x- 



the observed range for all of the objects in this sample. In 
contrast, the rest-frame 0.2 keV flux is derived from an ex- 
trapolation of the fitted data at all but the lowest redshifts, 
therefore the inferred flux at 0.2 keV depends upon the X- 
ray slope and Lo.2kcV may show an induced correlation with 
a x . 

To illustrate this problem, Z/2keV and I/o.2k c v are plotted 
as a function of a x in Fig 2. There is no apparent relation- 
ship between the Z/2kcV and the X-ray slope (q x ), but there 
is a strong correlation between a x and -Lo.2keV, i.e. when the 
X-ray spectrum is very soft, L .2kcv is high and when o? x is 
hard, Lo.2k c v is low. This is consistent with the anticipated 
selection effect. To verify that the trends in the observed 
data are similar to those in a random distribution, random 
samples of a x and L2keV were generated. The random data 
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were created by calculating the mean and standard devia- 
tions of the observed a x and L2k e v and using these to gen- 
erate random samples with a normal (Gaussian) probability 
distribution. The Lo.2k c v of this comparison sample is the 
sum of the random a x and Z/2kcV (recalling that Lo.2kcV and 
^2kcv are defined here as logarithmic quantities). 

The random data are compared with the observed quan- 
tities in Fig 2; both the Z/2k c v and Lo.2keV distributions with 
q x are similar in the observed and simulated datasets. This 
shows that the apparent correlation between q x and Lo.2k e v 
is consistent with the limitations of the sample's selection 
criteria and may not therefore reflect a 'true' property of 
AGN. 

The spurious dependence between Lo.2kcV and a x might 
also affect relationships with L .2kev and with a OB (the 
2500 A-to-0.2 keV ratio) and the significance of this effect 
is discussed where appropriate in Section 3.2. We do not 
expect any spurious dependences between the optical pa- 
rameters (i.e. a opt and L2500) and Z/2kcV or a x . 

2.2 Optical data 

The optical spectra were obtained over several observing 
runs with the Isaac Newton (INT) and William Herschel 
Telescopes (WHT) at the Observatorio del Roque de los 
Muchachos, La Palma. Two different instruments were used, 
the Faint Object Spectrograph (FOS) on the INT and the 
Intermediate-Dispersion Spectrograph and Imaging System 
(ISIS) on the WHT (these are indicated in Table 1). The 
FOS spectra cover a range of 3500A to 10000A with a res- 
olution of 15-20A FWHM in the red and 8-10A FWHM in 
the blue, while ISIS spectra cover 3000A to 9000A with a 
resolution of 3A FWHM in the red and 2A in the blue. All 
spectra were taken using a narrow slit and with the slit po- 
sitioned at the parallactic angle except where indicated in 
Table 1. 

CCD images of 23 of these AGN were also obtained 
at the Jacobus Kapteyn (JKT) and Nordic Optical (NOT) 
Telescopes and these were used to check for the amount of 
light typically lost around the narrow slit. For these AGN, 
the average ratio of fluxes measured from the CCD images 
to fluxes measured from the optical spectra was 1.23±0.12. 
Objects which had a strong galactic contribution in the CCD 
images were not included in this calculation. 

2.2.1 Optical slope 

Measurements of the optical continuum slope were made by 
fitting a single power-law to the spectrum, having first re- 
moved all emission and absorption features. The full wave- 
length coverage of the spectrum was used, but ignoring re- 
gions of poor data. No correction was made for any host 
galaxy contamination, although we expect that this effect is 
small in AGN with a redshift, z > 0.25. It has been demon- 
strated that the combination of Balmer continuum emission 
and blended Fell lines may produce a 'quasi'-continuum that 
may affect the measurement of the slope of any underlying 
continuum (e.g. Wills, Netzer and Wills 1985). In this analy- 
sis, no correction for the Balmer continuum or Fell emission 
has been made unless they have formed a distinct emission 
feature; in this case the features were removed before fitting 
the overall continuum. 



The best-fit observer-frame optical power-law indices 
(aiopt) are listed in Table 1. Indices were not measured from 
spectra not taken at the parallactic angle. Errors on the 
optical slopes are dominated by systematics and are difficult 
to determine for individual spectra. We have assessed the 
typical error expected on the slopes from the dispersion in 
the measured values for source observed more than once 
and estimate that the uncertainty in a op t is conservatively 
~ ±0.5. 

2.2.2 Luminosities 

The logarithm of the optical continuum luminosity in the 
quasar rest-frame at 2500A (L2500; also given in Table 1) 
was calculated from the best-fitting power-law model fit to 
the spectra. Luminosities of sources for which CCD data 
were also available were explicitly corrected for any light 
lost around the slit; these are indicated in Table 1. Where 
CCD images are not available, the luminosities have been 
increased by the mean factor of 1.23 (see Section 2.2). For 
spectra not taken at the parallactic angle, luminosities were 
calculated from a power law spectrum which was fitted to 
the red part of the spectrum only (these are also indicated in 
Table 1). Again, errors on L2500 are dominated by system- 
atics; taking into account possible variability, light losses 
around the slit and errors on the power-law fits, we esti- 
mate that uncertainties on optical luminosities are typically 
-50%. 

3 RESULTS 

In the remainder of this paper, we shall include a discussion 
of the optical/UV BBB and the soft X-ray excess in AGN 
continua, thus it is important to differentiate between these. 
Therefore we use the term 'BBB' to refer to the rise in the 
optical/UV continuum towards the blue, the 'soft X-ray ex- 
cess' to the excess above an extrapolation of the hard X-ray 
power-law at energies below ~1 keV, and the 'big bump' 
refers to an optical-to-soft X-ray component which incorpo- 
rates the two former features. 

We shall be describing much of the changes observed 
in the RIXOS optical and X-ray spectra in terms of the 
gradients of the continuum slopes. Different terminologies 
are commonly used in the literature for the slopes in the 
optical and the X-ray ranges, thus to avoid any confusion, 
we have adopted a single convention for both, i.e. using the 
terms "soft" and "hard". A "soft" slope falls towards high 
energies and has a relatively high energy index a (recalling 
here that the negative sign is implicit in our definition of 
a); in the optical region, a soft slope corresponds to a "red" 
continuum. A "hard" slope rises towards high energy and 
has a low or negative a; this corresponds to a "blue" slope 
in the optical. To describe changes in slope, we will use the 
terms "soften" (i.e. a increasing) and "harden" (where a is 
decreasing) . 

3.1 Sample properties 

3.1.1 X-ray power-law slope 

The average a x for the RIXOS AGN is 1.07±0.63 (la stan- 
dard deviation); this is harder than that of WF (1.50±0.48) 
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Table 1: Optical and X-ray continuum parameters. 
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Tabic 1 (continued): Optical and X-ray continuum parameters. 
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Tabic 1 (continued): Optical and X-ray continuum parameters. 
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97 79 +0.09 
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«opt 

(9) 

1.1 
0.2 
-0.1 
3.0 
1.0 
2.0 
0.5 
0.6 
0.4 
2.6 



(10) 

1.66 

1.40 

0.80 
-1.54 
-0.68 
-0.03 

1.78 

0.83 

1.20 

0.27 



Oios 
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1 0+°' 8 
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Oiox 

(12) 

1.40 
1.33 
1.32 
0.96 
1.18 
1.25 
1.38 
1.20 
1.03 
0.73 



+0.08 
0.13 
0.07 
-0.14 
0.09 
-0.13 
+0.09 
0.13 
+0.09 
0.13 
+0.05 
0.15 
+0.07 
0.14 
0.06 
0.15 
0.09 
0.13 
+0.05 
0.16 



(1) RIXOS field number (see Mason et al 1996); (2) RIXOS source number (Mason et al 1996) - the number in brackets is the 
radius of the extraction circle used for the X-ray data (in arcseconds) where it is less than 54 arcsec (see Section 2.1); (3) the 
telescope at which the spectrum was taken - see Section 2.2; (4) redshift; (5) Galactic column density (10 20 cm -2 ) - errors 
are ~10 per cent (see also Section 4.2.1); (6) log of the monochromatic optical luminosity at 2500 A (erg s _1 Hz -1 ) - error is 
estimated to be ~50 percent (Section 2.2.2); (7) log of the monochromatic X-ray luminosity at 0.2 keV (erg s Hz ) - errors 
are calculated from the 90% errors on the fits (Section 2.1.1); (8) log of the monochromatic X-ray luminosity at 2.0 keV (erg 
s~ 1 Hz~ 1 ) - errors are calculated from the 90% errors on the fits (Section 2.1.1); (9) energy index of the best-fitting power-law 
to the optical continuum - error is estimated to be ±0.5 (Section 2.2.1); (10) energy index of the best-fitting power-law to the 
X-ray data - errors are 90% (Section 2.1); (11) and (12) for definitions see Section 3.2.3 - errors are calculated from the quoted 
errors on L2500, ^o.2keV an< i ^2kev; 1 - CCD image taken; np - spectrum not taken at the parallactic angle. 
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Figure 3. The X-ray spectral index, a x , plotted as a function of 
redshift. 

and suggests that the WF sample has a bias towards AGN 
with strong soft X-ray excesses and thus perhaps also strong 
BBBs. The RIXOS mean a x is also harder than for the L94 
sample of optically-selected quasars (1.50±0.40). 

We find no significant dependence of q x on redshift 
for this sample of the RIXOS AGN (see Fig 3), neither is 
there any such dependence in the complete sample (Mit- 
taz et al. 1996). However the scatter in a x may lessen at 
high redshifts, e.g. for 2 < z < 3, the mean a x =0.85±0.09 
(standard deviation), which is similar to the mean slope of 



a x =0.89±0.25 measured at higher energies in low-z AGN 
by Comastri et al. (1992). This would be consistent with 
two competing trends, 1. a tendency for the soft X-ray ex- 
cess to be shifted out of the observed range as the redshift 
increases, and 2. a reduced sensitivity to absorption in the 
rest-frame at higher redshifts, as the most strongly affected 
regions (at energies below ~ 0.5keV) also move out of range. 



3.1.2 Optical continuum slope 

The best-fitting optical power-law slopes have been calcu- 
lated over the full range of the observed spectra and there- 
fore correspond to different rest-frame wavelength ranges at 
different redshifts. This must be taken into account when 
considering a op t, as previous studies have shown that the 
optical/UV continua of individual AGN soften towards the 
UV (e.g. Kinney et al. 1987; O'Brien, Gondhalekar & Wil- 
son 1988; Zheng & Malkan 1993), thus we might expect the 
RIXOS a op t distribution to show a general softening with 
redshift. The optical indices are plotted in Fig 4 as a func- 
tion of z, but at z > 0.25 (where galactic contamination is 
small), there are no changes with redshift within la. This 
appears to be contrary to the previous studies mentioned, 
but requires further investigation with spectrophotometric 
optical-to-UV coverage in the rest-frame to assess its signif- 
icance. 

The mean a opt for the RIXOS AGN is 0.92T0.76 (la- 
standard deviation). This is softer than the mean slope mea- 
sured for samples of optically-selected quasars, e.g. Neuge- 
bauer et al. (1987) found a median a opt of 0.2 while Francis 
et al. (1991) measured a median index of 0.3. This suggests 
that the optically-selected samples may be biased towards 
objects with relatively strong optical/UV BBBs. 
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Figure 4. The slope of the best-fitting power law to the optical 
continuum, a op t, plotted as a function of redshift. 



3.1.3 Ratio of optical to X-ray fluxes 

The ratios of optical (i.e. 2500 A) to X-ray flux at 2.0 keV and 
0.2 keV are parameterized using a ox (Tananbaum et al. 1979) 
and a s (Puchnarewicz et al. 1992) respectively. The a ox 
and a os are plotted as a function of redshift in Fig 5 (and 
listed in Table 1) and show no significant changes with z 
(although there is the suggestion of a possible increase in 
a os at z > 1.5). 

The mean a ox for the RIXOS AGN is 1.14±0.18 (la 
standard deviation). This is lower than for optically- 
selected objects, e.g. for the complete sample of PG quasars, 



the mean o? ox is 1.50 with a range of 1.0 to 2.0 (Wilkes 
et al. 1994). It is also lower than the mean which we calcu- 
late for the WF sample of X-ray-bright AGN, i.e. 1.37±0.23 
(la standard deviation). Note that the values of <a ox for 
WF have slightly different definitions to those used for the 
RIXOS objects because WF measure the optical continuum 
flux at 2675 A. For WF, we have calculated the optical-to-X- 
ray ratio using q ox (2675 A)=(Z/2675A--I'2kev)/2.635 (where 
i2675A is the logarithm of the rest-frame luminosity at 
2675 A). By calculating the a ox (2675 A) for the RIXOS 
AGN, we estimate that the mean difference between the 
two definitions of a ox [i.e. a ox -a ox (2675 A)] is +0.02. When 
taking into account the difference in the definitions, we find 
that the 'corrected' WF mean would be 1.39±0.23, increas- 
ing the discrepancy between the samples. 

The mean 2500 A-to-0.2 keV ratio (a os ) is 1.17±0.37 
(la); however, note that this value may be biased by the 
selection effect inherent in L .2kev (see Section 2.1.2). 

3.2 Correlations 

All correlations in this section have been evaluated using the 
Spearman rank-order correlation coefficient which is a non- 
parametric method insensitive to outlying points. The 'cor- 
relation probabilities' given in the following sections, fWr, 
are expressed as percentages and are equal to 100% minus 
the percentage probability of being drawn from a random 
sample. The correlation probabilities are given for all con- 
tinuum parameters in Table 2. 

One problem in searching for correlations in flux-limited 
samples of this kind, is that luminosities are dominated by 
the correction for the distance to the quasar, so that the 
continua of low-redshift objects have low luminosities while 
high-z objects tend to have high luminosities. This usually 
results in very strong correlations between luminosities in 
different bands and may not reflect a real dependence. In 
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Figure 5. The ratios of optical-to- X-ray luminosity at 2.0 keV and 0.2 keV, (a) o ox and (b) a s respectively, plotted as a function of 
redshift. Also shown are the mean indices in intervals of Az = 0.5 for < z < 2 and the mean for 2 < z < 3; vertical error bars represent 
the standard deviation (<r) on the means. 
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Table 2: Continuum correlation probabilities (p cor r; percent) 
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(-0.2 ) 52 
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Otoe 


(0.2 ) 


55 


(0.4 ) >99 


(-0.1 ) 2 
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(0.0 ) 4 


(-0.6 ) >99 




"ox 


(0.1) 
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(0.5 ) >99 


(0.3 ) 95 


(0.1) 3 


(-0.3 ) 97 


(0.4 ) >99 


(0.4 ) >99 



The correlation coefficients between two different luminosities were calculated after having adjusted for redshift (see Section 
4). 



order to allow for this, when comparing parameters which 
are both luminosities, each pair of parameters is divided by 
a normalizing factor which has the redshift dependence built 
into it. The factor used is the rest-frame luminosity equiva- 
lent to an observed flux of lx 1CP 15 erg cm~ 2 s _1 (arbitrarily 
chosen to match the observed-flux typically measured in the 
optical emission lines). In Table 2, correlation probabilities 
have been calculated after dividing the luminosities by this 
factor. 



3.2.1 Luminosity and X-ray slope 

We find no relationship between the X-ray spectral slope 
and the luminosity at 2 keV; the apparent correlation be- 
tween a x and L .2kcv is due to a selection effect (see Section 
2.1.2). There is no correlation between a x and the optical 
luminosity at 2500 A. 



3.2.2 Luminosity and optical slope 

Table 2 shows that the slope of the optical continuum is 
strongly anti-correlated with L2500, i.e. as the 2500 A lu- 
minosity increases, a opt becomes harder. This could be 
a consequence of measuring L2500 from optical power-law 
slopes which are measured in the observer's frame; the rest- 
frame flux at 2500 A does not fall within the observed wave- 
length range below 2=0.4 and must be extrapolated from 
the power-law fit at longer wavelengths. L2500 would then be 
high when a opt was hard which might influence the correla- 
tion. It may also be due to contamination by the host galaxy 
which would have the effect of softening a op t when the nu- 
clear emission (and therefore L2500) is low, indeed, Wilkes 
et al. (1994) found that starlight contamination was an im- 
portant factor for AGN with an L2500 lower than ~29. How- 
ever at 2=0.4, the rest-frame wavelength range over which 
Q op t is calculated lies at ~2500-6400 A where galactic con- 
tamination is already relatively low and should have only a 
small effect (if any) on the measurement of the optical slope. 

Thus to eliminate these effects, we have excluded the 
2 < 0.4 objects and find that while the anti-correlation does 
disappear for z > 1.0 (rs=-0.1; p C orr=2%), it remains strong 
at 0.4 < 2 < 1.0 (r s =-0.6 and p corr =99.98%; see also Fig 6). 
A similar effect is seen in the correlations between a opt and 
^o.2kcV and between a opt and L 2 k c v (i.e. the data are cor- 
related for 0.4 < 2 < 1.0, but not at z > 1.0), although the 
latter dependences may be induced by strong correlations 
between L2500, £<o.2keV and L 2 k c v (see Table 2). 
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Figure 6. The optical spectral index, o pt , plotted as a function 
of optical luminosity at 2500 A, 1/2500- The AGN arc plotted 
according to their redshift, i.e. AGN with z < 0.4 are plotted 
as open squares, 0.4 < z < 1.0 arc plotted as solid squares and 
z >1.0 as asterisks. 

This might be a luminosity, rather than a wavelength- 
dependent effect, but it is difficult to distinguish between 
these due to the strong redshift dependence in the luminosi- 
ties. Nonetheless, the correlations imply that at relatively 
low bolometric luminosity, the A > 5000 A region hard- 
ens as the luminosity increases, while at higher bolomet- 
ric luminosities, the continuum slope at shorter wavelengths 
(A < 5000 A) is independent of luminosity. A similar effect 
was reported by Zheng & Malkan (1993), who found that the 
rest-frame 6500-3100 A continua of their quasars and Seyfert 
Is also harden significantly with luminosity (Aa=l.l) while 
at shorter wavelengths (3100-1250 A in the rest- frame), the 
change in slope is smaller over the same luminosity range 
(Aa=0.5). To explain this, they suggest that either (1) the 
BBB shifts from the UV into the optical as the luminosity in- 
creases, or (2) that the BBBs are stronger in high-luminosity 
AGN. 

3.2.3 Optical-to-X-ray luminosity ratios 

We find evidence for a strong correlation between a ox and 
Q x (p C orr=99.8%; see Fig 7). The 2500 A luminosity is mea- 
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sured independently of the X-ray parameters and L2keV is 
not correlated with a x , thus we conclude that this is not 
due to any selection effects. Similar correlations have also 
been reported by WF and L94. This relationship suggests 
that when the 2500 A-to-2 keV ratio is high, the 0.2-2 keV 
slope is relatively soft so that it rises steeply towards the 
EUV, and the overall 2500 A-to-2 keV continuum is convex. 
Conversely, when a ox is low, a x is hard (i.e. it falls towards 
the EUV) and the 2500 A-to-2 keV continuum is concave. 
Assuming that a ox is an indicator of the strength of the 
BBB (because the BBB is strong at 2500 A but negligible 
at 2 keV) and that ev x measures the strength of the soft X- 
ray excess, then this implies that the optical/UV BBB and 
the soft excess may be part of the same component. 

It is important to be mindful of possible selection effects 
when considering such correlations, particularly in the X-ray 
data where the spectral resolution is low, and the response 
of the detector favours particular energies within the nom- 
inal bandpass. For example, Ulrich & Molendi (1996) have 
cautioned that the 2 keV flux inferred from ROSAT spectra 
is coupled to the spectral slope. They suggest that this may 
induce the correlations seen by WF and L94 because 2 keV, 
the reference energy for q ox , is far removed from the cen- 
troid energy of the counts detected by the ROSAT PSPC 
for a typical low-redshift Seyfert, which is ~0.4 keV for a 
slope, q x ~1.6. However, there are a number of reasons why 
this effect is unlikely to induce the correlations seen in the 
RIXOS AGN. Firstly, the typical error bars, which reflect 
the coupling between spectral normalization and slope in 
the size of the uncertainty in q ox , are significantly smaller 
than the observed range of the correlation. Secondly, the 
X-ray spectra of objects in our sample are typically harder 
than those in WF and L94 and cover a wider range in a x ; 
both of these factors will reduce the Ulrich & Molendi effect. 
Thirdly, the RIXOS AGN extend to much higher redshifts 
than WF and L94 moving the equivalent of the 2 keV flux 
in the quasar's rest-frame closer to and below the centroid 
energy of the photons registered in the detector frame, yet 
the correlation is still seen in the sample of RIXOS objects 
above 2=1 (see Fig 7). 

We also find evidence of a correlation between L . 2 k e v 
and a ox (j5corr=95%), i.e. when the optical-to- X-ray ratio is 
high, the luminosity at 0.2 keV is also high, although this 
may be an induced correlation as both L .2kcv and a ox are 
also correlated with q x (and the £o.2keV-c*x correlation is 
due to a selection effect; see Section 2.1.2). 

Table 2 also indicates a correlation between L2500 and 
a ox , but this is expected since the L 2 k e v distribution is ran- 
dom with q x (Section 4.2.1) and L2500 is closely related to 
q ox (i.e. 1/2500=2. 605a ox +£2kcv)- However, significant de- 
pendences have been reported elsewhere between luminosity 
and a ox (e.g. Reichert et al. 1982; Avni & Tananbaum 1986; 
Wilkes et al. 1994). The correlation between L2500 and a os 
is consistent with the selection effect discussed in Section 
2.1.2. 

3.2.4 Optical and X-ray slopes 

There is an anti-correlation between a opt and a ox , imply- 
ing that when the 2500 A-to-2 keV ratio is high, the optical 
slope is relatively hard. There is also a weak (p C orr=68 per 
cent) anti-correlation between the optical and X-ray slopes 
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Figure 7. The optical-to- X-ray ratio, oj ox > plotted as a function 
of X-ray spectral index, q x . The data are plotted according to 
the redshift of the AGN, i.e. z < 0.5 are plotted as diamonds, 
0.5 < z < 1.0 are plotted as filled boxes and z > 1.0 as crosses. 







oz<0.5 







. 0.5 <z< 1.0 : 






Hz > 1.0 




EB 













■ 








°o ' o 


: 

<a ~ 




a IN 


: 




m ^% ' 




% ; 




EST E 


f^O - 




■ i^^^^fj 






EE H|f| 












EB ■ 






Median 90% error 


1 1 . .- 


-2 


-1 1 


2 3 



Figure 8. Optical power-law slope, « op t, plotted as a function 
of the X-ray power-law slope, a x . The data are plotted according 
to the redshift of the AGN, i.e. z < 0.5 are plotted as diamonds, 
0.5 < z < 1.0 are plotted as filled boxes and 2 > 1.0 as crosses 
within boxes. 

when the sample is considered as a whole, however there 
appears to be a redshift dependence in the relationship (see 
Fig 8). At low redshifts, a opt is likely to be affected by the 
host galaxy and thus may be soft relative to higher-z AGN 
(cf. Fig 4); indeed Fig 8 demonstrates that the z < 0.5 
AGN are systematically displaced to softer values of a opt 
relative to the higher redshift objects. When the data are di- 
vided into redshift bands, the correlation strengthens, e.g. at 
z < 0.5, p C orr=80 per cent (rs=-0.4) while for z > 0.5, p CO rr 
also rises to 80 per cent (rs=-0.3). However, the strongest 
correlation is found in the high-z AGN, i.e. for z > 1.0, 
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Pcorr=96 per cent (rg=-0.5). The degree of correlation (as 
measured using the Spearman rank-order coefficient) is af- 
fected by the tendency of the data to cluster around a x =l 
(see Fig 3) and the relative paucity of very hard and soft 
q x objects which is an expected consequence of the RIXOS 
selection criteria. Samples with a more homogeneous dis- 
persion in a x and improved measurements of the intrinsic 
slopes will provide a more definitive answer. 

Nonetheless, even given the relatively basic methods for 
measuring q x and a op t, the data for these AGN show com- 
pelling evidence for a correlation between the two slopes. 
The correlation implies that when the optical continuum is 
hard, the soft X-ray slope is soft (and vice-versa) so that 
both the optical and X-ray slopes fall and rise towards the 
EUV together. 

This is the first time that a relationship between the 
optical and X-ray slopes has been found. L94 found no 
correlation in their PG quasar sample, but since the RIXOS 
data suggest that a op t changes only very slowly with a x (see 
Fig 8), a sample covering a wide range in a x and/or a opt 
is essential to identify any trend; the L94 AGN were UV- 
excess quasars whose range in a opt (-0.2 to 0.9) may have 
been too small for this correlation to be detected. 

Together with the dependence of a ox on a x , these corre- 
lations are additional evidence linking the optical/UV BBB 
to the soft X-ray excess. In short, when a x is soft, the opti- 
cal continuum (a op t) is hard and the optical-to- X-ray ratio 
(i.e. q ox ) is high. Conversely, when a x is hard, a opt is soft 
and a ox is low. 

3.3 Mean optical-to-X-ray continua 

To illustrate the changes in the shape of the optical-to-X- 
ray spectrum, we have calculated three 'mean' spectra at 
different values of a x , covering 6000 A to 2 keV. The mean 
a ox , Q os and a opt were calculated in three bins centred on 
a x =0.25, 0.95 and 1.65 and each covers an interval Aa x =0.7. 
Only AGN with redshifts in the range 0.25 < z < 3.25 
were used to avoid possible host galaxy contamination in the 
optical at low redshifts and the highest-z object whose X-ray 
and optical continuum slopes have relatively high errors. 

The mean indices and the errors on the means are shown 
in Fig 9a-c and listed in Table 3. These were then used 
to construct the mean 'spectra' at o? x =0.25, 0.95 and 1.65 
(Fig 9d; they have been normalized at 2.0 keV) and illustrate 
more clearly the changes in the spectra suggested by the 
data. When q x is very soft, a op t is hard and the spectrum 
is "convex" (indicated by the dashed line in Fig 9d), while 
when a x is hard, a op t is soft and the overall spectrum is 
"concave" (the dot-dashed line in Fig 9d). 



4 MODELLING THE CONTINUUM 

Recent work on the continua of Seyfert Is and quasars has 
contributed much to our understanding of AGN spectra, es- 
pecially with regard to the possibility of a single 'big bump'. 
Both WF and L94 identified a strong o? ox vs <a x correlation; 
WF suggested that changes in the flux of a big bump rel- 
ative to an invariant underlying component could explain 
the spectral differences from source to source and that the 
shape of the big bump from the UV to soft X-rays showed 
little change. 



Table 3. Mean optical and X-ray spectral indices 
(0.25 < z < 3.25) 



-0.1<a x <0.6 

a ox 1.10±0.05 
a os 1.62±0.12 
a op t 1.50±0.33 



0.6<a x <1.3 

1.13±0.02 
1.21±0.04 
0.72±0.10 



1.3<a x <2.0 

1.21±0.02 
0.99±0.04 
0.72±0.10 



Errors quoted are the errors on the mean (i.e. ct/(N-1), 
where N is the number of data points). 

The overall changes in optical-to-X-ray continua of the 
RIXOS AGN (illustrated in Fig 9d) could also be inter- 
preted as being due to the sum of two components; a big 
bump spectrum, whose flux varies significantly from source 
to source, and a relatively invariant underlying continuum 
(i.e. a similar model to that proposed by WF). If this model 
is to satisfy the full range of the observed RIXOS spectra, 
the underlying continuum must be concave with a break be- 
tween 2500 A and 0.2 keV (see the mean spectrum indicated 
by the dot-dashed line in Fig 9d; we stress that we have no 
information about the shape of spectrum between 2500 A 
and 0.2 keV). 

However, an underlying continuum of this form is dif- 
ficult to reconcile with observations at longer and shorter 
wavelengths. The slope of the underlying continuum in 
soft X-rays (i.e. based on the RIXOS AGN with the hard- 
est a x ) would be equivalent to q x ~0, which is significantly 
harder than the indices of Seyferts and quasars at higher 
energies [~2-50 keV; for Seyferts a x =0.7 while for quasars, 
q x is closer to 1.0 with small dispersions about the mean 
(Mushotzky, Done and Pounds 1993 and references therein)]. 
Also, an extrapolation of the rather soft optical continuum 
to longer wavelengths (equivalently a op t~2 when a x ~0) 
predicts a very strong flux in the IR, whereas typically 
the continua of AGN spectra are generally harder with an 
a ~1 extending down to ~100/mi (for examples see Elvis 
et al. 1994). 

It is possible to have an underlying continuum of this 
shape and remain consistent with the IR and hard X-ray 
continua by assuming that the underlying continuum is a 
simple power-law which is absorbed predominantly in the 
optical/UV and soft X-ray regions, e.g. by cold gas and 
dust. However, it is in practice difficult to envisage a phys- 
ical model where the nuclear power-law source lies behind 
significant amounts of cold, dusty gas, leaving the big bump- 
emitting region bare. 

Thus we conclude that the changes in continuum shape 
for the RIXOS AGN cannot simply be due to either 1. a 
continuum which is the sum of a power law and a big bump 
with different relative strengths, and which is unabsorbed; 
or 2. a continuum which is the sum of an absorbed power 
law and an unabsorbed big bump. 

Neither are the observed correlations consistent with 
a big bump component which 'shifts' between the optical 
and soft X-ray regions. In this case, when the big bump 
spectrum lies at relatively low frequencies (i.e. it is strong 
in the optical but weak in soft X-rays), we would expect 
a hard a op t and a hard q x ; conversely when the spectrum 
lies at relatively high frequencies (i.e. weak in the optical 
but strong in soft X-rays), both a opt and a x would be soft. 
Thus we would expect to see a correlation of a op t with a x 
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Figure 9 The calculation of three 'mean' RIXOS AGN spectra from 6000 A-to-2 kcV. (a)-(c) Data in the range z=0.25 to z=3.25. Also 
plotted as horizontal bars are the means in each of three equally spaced o x bins for these data (vertical bars indicate the errors on the 
means; see Table 3). (a) a ox plotted against ct x ; (b) a OB plotted against ct x ; (c) a op t plotted against ct x . (d) Three mean continuum 
spectra calculated from the means shown in (a)-(c). Only the separate optical and X-ray slopes and their relative normalization are 
known; we have no information in the 2500 A-to-0.2 keV region. The spectra may be identified by the broken lines joining the optical 
and X-ray regions; i.e. a x =0.25 is indicated by the dot-dash line (also referred to as the "concave" spectrum), o x =0.95 is indicated by 
the dotted line and a x =1.65 (referred to as the "convex" spectrum) is shown by the dashed line. 



rather than the anti-correlation observed. Examples of shift- 
ing spectra are AD models with different black hole masses 
and bremsstrahlung components with a range of tempera- 
tures. 

4.1 Cold absorption 

Instead we propose that the changes in continuum shape 
which are responsible for producing the observed correla- 
tions, may be primarily due to absorption of the entire in- 
trinsic spectrum by cold gas and dust. The cold gas would 
be responsible for absorption in soft X-rays while the dust 
would be the source of absorption in the optical and UV. 



The intrinsic spectrum may be the sum of a power-law and 
a big bump. The big bump temperature, power-law slope 
and big bump-to-power-law continuum flux ratio may also 
vary from source to source, introducing scatter into the data, 
however it is the absorption which results in the observed 
correlations. 

To test whether this crude model is consistent with the 
observed spectra, the optical and X-ray spectral slopes and 
ratios were simulated by constructing a model for the intrin- 
sic spectrum and modifying this for the effects of absorption 
by cold gas and dust. Throughout the remainder of this sec- 
tion, we investigate whether this model can predict the data 
which is actually observed, while the physical consequences 



Properties of RIXOS AGN: I - continua 13 



of such a model for AGN are discussed in Section 5. This 
is intended to be a preliminary assessment and consistency 
check of the model; a more rigorous treatment awaits a fu- 
ture paper. 

The intrinsic spectrum in the model was assumed to be 
the sum of an underlying continuum and a big bump and 
its flux was fixed while the amount of absorption was var- 
ied. A simple power law with a slope of a=1.0 was adopted 
for the underlying optical-to-X-ray continuum; the chosen 
model for the big bump is discussed in Section 4.1.1. In 
soft X-rays, the Morrison and M c Cammon (1983) absorp- 
tion cross-sections were assumed while in the optical and 
UV, the reddening curve of Cardelli, Clayton & Mathis 
(1989) was used. A Galactic dust-to-gas ratio [where E(B- 
V) = 1.0 corresponds to an Ah of 6xl0 21 cm~ 2 ; Ryter, Ce- 
sarsky & Audouze 1975; Gorenstein 1975] was assumed for 
the dust abundance. The model a opt was calculated by fit- 
ting a power-law to the model spectrum; the 2200 A dust 
absorption feature was ignored in the fitting as this is rarely 
observed in the UV spectra of AGN (Cardelli & Clayton 
1991; Mathis & Cardelli 1992). The X-ray spectral index 
was derived by first binning the model spectrum into three 
'colours' which reflected the binning used in the analysis of 
the data (see Section 2.2), then fitting a power-law to these 
colours. 

Because the AGN spectra sample such a wide range in 
z, it is necessary to include the effects of redshift on the 
model. The optical and X-ray slopes have been measured in 
fixed ranges in the frame of the observer, therefore they cor- 
respond to different parts of the intrinsic quasar spectrum, 
depending on z, and this must be taken into account when 
trying to reproduce the data. Therefore the models have 
been calculated at different redshifts (z—0.2, 0.5, 1.0, 1.5 
and 2.5) to check that the predictions are consistent with 
the data. Both the optical and X-ray ranges were corrected 
by a factor of 1 + z before the model indices were calculated. 

4-1.1 Choosing a model for the big bump 

To find an appropriate model for the big bump, we have 
compared the "convex" spectrum in Fig 9d (i.e. the mean 
a x =1.65 RIXOS spectrum, which we assume to be domi- 
nated by the big bump) with several different models in- 
cluding ADs and bremsstrahlung spectra. The convex spec- 
trum is much broader in wavelength than for many types of 
AD model, whether geometrically-thick (e.g. Madau 1988) 
or thin (e.g. Sun and Malkan 1989), although it does com- 
pare well with the modified Czerny & Elvis (1987) model of 
an AD which is surrounded by a hot (T = 10 8 K) corona. 

A simple thermal bremsstrahlung spectrum with a tem- 
perature, Tbrom=l-2x 10 K also compares well with the con- 
vex spectrum; bremsstrahlung models have previously been 
used to fit the big bumps of other AGN with soft excesses, 
e.g. PG1211+143 was fitted with a T brcm =8xl0 5 K spec- 
trum (Barvainis 1993) and E1346+266 compared well with a 
bremsstrahlung model with Tb rcm =7xl0 6 K (Puchnarewicz 
et al. 1994). 

The AD-plus-corona and bremsstrahlung models are 
plotted in Fig 10 with the mean o? x =1.65 'observed' spec- 
trum (an q=1.0 power-law has been added to the model 
spectra). In the subsequent modelling, for simplicity, we 
have used a thermal bremsstrahlung spectrum for the big 
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Figure 10. A comparison of the mean convex RIXOS spectrum 
(i.e. the mean a x =1.65 spectrum from Fig 9d; plotted here as a 
solid/long-dashed line ) with two models for the BBB and soft 
X-ray excess; 1. the Czerny & Elvis (1987) model of an accretion 
disc surrounded by a hot (T = 10 s K) corona (dotted line), and 
2. a Tb rcm =1.2x 10 6 K simple bremsstrahlung spectrum (dashed 
line). An a x =1.0 power-law has also been added to both models. 

bump, although the AD-plus-corona provided an equally ac- 
ceptable alternative. 

4.2 Finding the best fit 

In this section, we assess the ability of the model to re- 
produce the observed correlations between q x and a ox and 
between a x and a op t ■ While we have searched for the best- 
fitting models over a wide range of parameter space (i.e. in 
Tbrcm and the relative normalization of the bremsstrahlung 
component), the model has not been allowed to converge 
formally, therefore the presence of other minima in the data 
are possible. This approach was considered adequate for the 
purpose of a first evaluation of the model. 

In searching for the model which provided the best rep- 
resentation of the data at each redshift, we changed two of 
the parameters, Tb rcm and the flux in the bremsstrahlung 
component relative to the power law. The model was calcu- 
lated over a grid of Tb rem and big bump-to-power-law flux 
ratios and at each point in the grid, T brem and the flux ra- 
tio were fixed while the intrinsic A^h was varied. The a ox 
versus a x and a op t versus a x distributions were fitted simul- 
taneously and the goodness-of-fit was evaluated using the xt 
statistic (taking into account the errors on a ox and a op t; see 
Table 1). The 'best-fits' described throughout the remainder 
of the paper represent the points in the grid with the lowest 
xt- The a os versus a x distribution is biased by a selection 
effect and may not be a real physical dependence (Section 
2.1.2), therefore it was not included in the fits (although the 
predicted model has been compared with the data in the 
plots as a check for consistency). The models were only fit- 
ted to the data appropriate to the redshift assumed, i.e. the 
models at 2=0.2, 0.5, 1.0, 1.5 and 2.5 were fitted to the data 
in the 2=0-0.25, 0.25-0.75, 0.75-1.25, 1.25-1.75 and 1.75-3.5 
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Figure 11. A comparison of the RIXOS AGN data with the model which assumes the sum of a bremsstrahlung spectrum and an a x =1.0 
power-law and which is absorbed by cold gas and dust. The flux in the power-law and the bremsstrahlung spectrum remain constant 
while the absorbing column density is varied; (a) illustrates the range of continua predicted by the 2=1.0 model: the thick line is the 
unabsorbed spectrum while each curve plotted as a thin line represents absorption by an additional 6xl0 20 cm~ 2 . In (b)-(d), models are 
plotted for z=0.2 (dot-dashed line), 2=0.5 (dashed line), z=1.0 (solid line) and 2=1.5 (dotted line). Model spectral indices are compared 
with the data in (b) where c? ox is plotted against a x ; (c) a OB is plotted against a x , and (d) a op t is plotted against c? x . In all cases, the 
observed data are plotted as diamonds. 



ranges respectively. 

The 'best-fitting' models at all redshifts except 2=2.5 
are plotted in Fig 11 and show a good agreement with the 
data, considering the simplistic assumptions made. Parame- 
ters of the models are listed in Table 4; the normalization of 
the bremsstrahlung component relative to the power-law is 
expressed as Ro.2keV, the ratio of flux in the bremsstralung 
spectrum to flux in the power law at 0.2 keV. 

There are very few differences between the best-fitting 
a ox vs a x models as the redshift changes, which is consistent 
with the observations (see Figs 4, 5 and 7, and Section 3.1). 
The model does predict however that for any given a x , a op t 
should be relatively soft at higher-z, whereas the observed 



values for a op t tend to be softer at low-z (i.e. z < 0.5; al- 
though this may be due to contamination by the host galaxy 
for which no corrections have been made) . The slopes of the 
predicted a OB versus a x models are harder than the observed 
data; this may be a consequence of the selection effect dis- 
cussed in Section 2.1.2. 

At higher redshifts (z > 1.75), the range of a x is rela- 
tively small and the relationships between the indices may 
not be the same (Section 3.1), therefore we discuss these 
objects separately. When z is high, the model predicts that 
the effects of absorption in soft X-rays are lessened (because 
the strongly absorbed part of the spectrum has been shifted 
out of the observed range), while in the optical, the absorp- 
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Table 4 Parameters of the best-fitting models to the 
RIXOS AGN continua (plotted in Figs. 11 and 12) 



z 


-Ma rem 


Ro.2keV 


N H 




(i) 


(2) 


(3) 


0.2 


1.7xl0 6 K 


10 


1.6 


0.5 


1.7xl0 6 K 


20 


4.0 


1.0 


1.3xl0 6 K 


20 


6.0 


1.5 


1.2xl0 6 K 


30 


7.0 


2.5 


1.0xl0 6 K 


40 


4.0 



(1) Temperature of the bremsstrahlung component. (2) 
Ratio of flux in the bremsstrahlung component to flux 
in the power-law component at 0.2 keV. (3) Maximum 
column density of the gas used to produce the model in 
units of 10 21 cm~ 2 . 

tion has a greater effect (because the heavily-absorbed UV 
is shifted into range). Therefore, the scatter in a x at high- z 
would be reduced while the range in a opt and a ox would 
remain relatively wide. In Fig 12 we compare the a ox ver- 
sus cv x and a opt versus a x distributions with the best-fitting 
model at z—2.5 (parameters of the model are listed in Tabic 
4) and we find that the high-redshift data are consistent with 
the predictions of the model. Thus the model can predict 
the observed data from z—Q to ~3. 

4-2. 1 Errors on the Galactic Nu 

The measurements of a x used in this analysis have already 
been corrected for A^Gai, therefore any absorption implied 
by the modelling is additional to this. In principle, this 
could be caused by errors on the assumed A^HGai rather 
than absorption external to our Galaxy. To test whether 
the changes in the continuum slopes of the RIXOS AGN 
could be due to errors on A^Gai, we have compared the A^h 
distribution implied by the model with the typical uncer- 
tainties on the Stark et al. (1992) measurements of NhgbI- 



By default the model assumes that the gas and dust absorp- 
tion occurs in the rest-frame of the quasar; to measure the 
equivalent Galactic column, the model was modified to as- 
sume that the absorption is entirely within the frame of the 
observer. 

We found that the typical additional observer-frame col- 
umn density suggested by the model would be ~ 2 x 10 20 
cm~ 2 (and this model also provided a poorer representation 
of the data, especially at high-z). Elvis, Lockman & Wilkes 
(1989) calculated a la error on the Stark et al. measure- 
ment of AfHGai of 11 per cent for pointings at high Galactic 
latitude, implying typical errors on A^Gai for the RIXOS 
AGN of ~ 2 x 10 19 cm~ 2 ; this is significantly lower than the 
additional 'Galactic' columns implied by the model. Thus 
we conclude that the variations in continuum shape are not 
simply a reflection of the errors on the measurements of 

A^HGal- 

4.3 Range of big bump temperature and strength 

Although the correlations between a ox and q x and between 
Q op t and o? x are significant, there is still a large degree of 
scatter in the data. While the modelling suggests that the 
correlations may be due to absorption, the scatter may re- 
flect a range in other parameters such as the temperature 
and relative normalization of the big bump. Therefore, by 
calculating models for different values of Tb re m and the big 
bump-to-power-law ratio and comparing these with the ob- 
served distributions, we can find the approximate ranges of 
temperature and normalization implied by the model. 

4-3.1 Tbrc m 

To investigate the range of bremsstrahlung temperature, dif- 
ferent models were calculated at low and high Tb rC m until 
they fell outside the regions occupied by the distributions 
of the spectral indices (i.e. a ox versus a x and a opt versus 
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a x ). The normalizations and column densities were adjusted 
where necessary to provide the best comparison with the 
data. 

The limits on the bremsstrahlung temperature implied 
by the model at a redshift of 1.0 are shown in Fig 13 and 
listed in Table 5, and show that a wide range in temperature 
can be accommodated (from 10 s to 3xl0 6 K). At lower tem- 
peratures, the bremsstrahlung spectrum does not reach far 
enough into soft X-rays to reproduce the soft excess, while 
at higher temperatures, the a ox is too low for any given a x , 
as the big bump component shifts too far out of the optical 
range. Both upper and lower limits are constrained by the 
Qox versus a x distribution; the scatter in the observed a opt 
versus a x data falls outside these limits (at low and high 
Q op t) but because of the relatively high errors on the mea- 
surements of Qf op t we are unable to determine whether this 
is physically significant. 

4-3.2 Ratio of big bump to power-law flux 

A similar procedure was followed to estimate the range of 
the big bump-to-power-law flux ratio implied by the model. 
In this case, r brC m was fixed while the normalization of the 
bremsstrahlung spectrum was varied until the model fell out- 
side the regions defined by the data. Only the column den- 
sity was adjusted where necessary to extend the model to 
hard a x . 

The results (at 2 = 1) are shown in Fig 14 and are 
compared with the best-fitting model (parameters are listed 
in Table 6) . When the bremsstrahlung normalization is zero, 
there is no soft X-ray excess or BBB and only sources with 
a relatively hard a x and soft a op t can be predicted by the 
model. At the other extreme when Ro.2keV is very high, 
there is little change in the model for the predicted a ox and 
Q os versus a x distributions, but at any given a op t, the X- 
ray slope is much softer. This is because the bremsstrahlung 
dominates throughout the optical/UV so that the slope is 
unchanged by the flux in the big bump (except for very low 
normalizations), whereas in soft X-rays, the bremsstrahlung 
spectrum only dominates at relatively low X-ray energies 
(i.e. less than ~0.5 keV) and therefore its strength relative 
to the power-law has a much greater effect on the measured 
X-ray slope. 

The model thus suggests that when Tb rcm =1.25x 10 6 K 
(at z—1.0), Ro.2keV can vary from to ~1000. The upper 
limit is constrained by the a opt versus a x distribution (al- 
though we note that the absorption required is correspond- 
ingly high to predict any given a x ) but there is no lower 
limit (i.e. the model can accommodate AGN without big 
bumps) . 

5 DISCUSSION 

We have presented the optical and X-ray continuum parame- 
ters of 108 X-ray selected AGN from the RIXOS survey (Ma- 
son et al. 1996). From the correlation analysis of this sam- 
ple, we find significant evidence for relationships between 
q ox and a x and between a op t and a x (especially at z > 1.0). 
These suggest that as the 0.2-2 keV spectrum softens, the 
2500 A-to-2 keV ratio increases and the slope of the optical 
spectrum hardens. The trends apply over a wide range of a x 
for soft and hard X-ray slopes (i.e. from q x ~0 to 2.5) and 



Table 5 Range of bremsstrahlung temperature at 
2=1: parameters of models shown in Fig. 13 



^brcm 


Ro.2keV 


JVh 


(i) 


(2) 


(3) 


lxl0 5 K 


10* 


3.0 


1.3xl0 6 K 


20 


4.0 


3xl0 6 K 


30 


5.0 



(1) Temperature of the bremsstrahlung component. 

(2) Ratio of flux in the bremsstrahlung component 
to flux in the power-law component at 0.2 keV. (3) 
Maximum column density of the gas used to pro- 
duce the model in units of 10 21 cm -2 . All mod- 
els are calculated at 2=1.0. ' Ratio calculated at 
0.02 keV; the 0.2 keV ratio (R . 2 kev) is lxl0~ 7 . 

Table 6 Range of bremsstrahlung normalization at 
2=1: parameters of models shown in Fig 14 



-^brcm 


Ro.2keV 


iV H 


(i) 


(2) 


(3) 


1.3xl0 5 K 


1000 


10.0 


1.3xl0 6 K 


20 


5.0 


1.3xl0 6 K 





4.0 



(1) Temperature of the bremsstrahlung component. 

(2) Ratio of flux in the bremsstrahlung component 
to flux in the power-law component at 0.2 keV. (3) 
Maximum column density of the gas used to produce 
the model in units of 10 21 cm~ 2 . All models are 
calculated at 2=1.0. 

suggest that the overall 6000 A-to-2 keV spectrum changes 
from convex to concave as o? x hardens. 

WF and L94 have also identified a ox versus a x corre- 
lations in their relatively soft X-ray samples and WF sug- 
gested that it is due to changes in the strength of a big 
bump relative to an underlying continuum. However, this 
explanation is not appropriate for the RIXOS sample (which 
extends to harder X-ray slopes than WF) because it implies 
that the underlying continuum from 6000 A to 2 keV must be 
concave (see e.g. the dot-dashed line in Fig 9d) and a spec- 
trum of this shape is difficult to reconcile with observations 
at longer and shorter wavelengths (Section 4) . 

We have illustrated the main differences between the 
WF, L94 and RIXOS AGN samples, by plotting a ox as a 
function of q x in Fig 15 (these have been corrected for the 
differences in the definition of a ox for the three samples, 
i.e. the WF and L94 a ox s have been increased by 0.02 and 
0.04 respectively; see Section 3.1.3). The figure shows that 
the WF objects are clustered towards high <a ox and high <a x 
relative to the RIXOS AGN; a similar and more extreme 
pattern is seen for the L94 objects. This would have led 
WF and L94 to draw different conclusions regarding the na- 
ture of the optical-to- X-ray continua in AGN, e.g. the WF 
data do not require the effects of absorption perhaps be- 
cause they observed very few highly absorbed sources (with 
q x <1), while L94 suggested that sources with a soft a x were 
weak in X-rays (i.e. a ox was high when a x was soft), yet Fig 
15 shows that in both the WF and RIXOS samples, there 
are also many sources with a much lower a ox when a x is soft 
(i.e. the 2 keV flux is high relative to the optical and the X- 
ray spectrum is soft), implying a true soft X-ray 'excess' in 
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Figure 13. The range of 7b rem at a redshift of 1.0 which is implied by the scatter in the observed data. Three different models are 
shown: for the dot-dashed line, 7t, rem =10 5 K; for the solid line, Tj-, rcm =1.25xl0 6 K; for the dot-dot-dot-dashcd line, Tt, rem =3x 10 6 K. 
All three models are compared with (a) the observed ct ox versus ct x distribution; (b) the observed o os versus ct x distribution; (c) the 
observed a op t versus o x distribution. The three unabsorbed model spectra are plotted in (d). 



these AGN. 

These differences are largely a consequence of the selec- 
tion methods, i.e. the PG quasars used by L94 were selected 
by the strength of their UV excess, thus they tend to be 
optically-bright and -hard, while the WF objects are bright 
in soft X-rays (0.08-2.4 keV) and therefore X-ray-bright and 
generally X-ray-soft. The RIXOS sample is itself restricted 
by its selection criteria (i.e. by the strength of the 0.4-2.0 keV 
emission), therefore its properties can only be said to be 
typical of objects which are relatively strong in medium X- 
rays (i.e. at ~2 keV in the quasar rest-frame). Ultra-soft X- 
ray objects like E1346+266 and RE J1034+396, which have 
very soft ROSAT spectra yet little or no flux above ~1 keV 
(Puchnarewicz et al. 1994; Puchnarewicz et al. 1995), would 
be selected against in this sample for instance; as would 
AGN which are heavily absorbed through cold gas. 



As an explanation for the observed trends in the RIXOS 
AGN, we have proposed that the correlations are due to the 
effects of moderate absorption by cold gas and dust. It is 
possible that the results are also reproducible by a dusty 
warm absorber (see e.g. Mathur et al. 1994; Brandt, Fabian 
and Pounds 1996), but the additional parameters required 
for the fits are not justified by the quality of these data. 
This hypothesis has been tested using simple combinations 
of power laws and bremsstrahlung spectra which have been 
modified for the effects of dust and gas absorption. We find 
that the model is consistent with the data up to z ~ 3 (this 
is effectively the limit of the sample). 

However, the model must also be judged by its abil- 
ity to reproduce the spectra of AGN which lie outside the 
RIXOS parameter space, such as ultra-soft X-ray sources 
and the PG quasars (see Fig 15), which are relatively quiet 
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Figure 14. An illustration of the effect of changing the strength of the bremsstrahlung component relative to the underlying power law. 
Three different models (with z=1.0) are shown and compared with the data: for the dot-dot-dot-dashed line, there is no bremsstrahlung 
component; for the solid line, the -f?o.2fceV=20; for the dot-dashed line, -Ro.2fceV=1000. AH three models have Tb rem =1.25xl0 6 K and 
are compared with (a) the observed ct x versus ct x distribution; (b) the observed a OB versus a x distribution; (c) the observed a op t versus 
o x distribution. The three unabsorbed model spectra arc plotted in (d). 



in X-rays. We have fitted our model to the L94 sample of 
PG quasars, finding a best-fit when T brcm =6xl0 5 K and 
Ro.2keV=W, which reproduces the optical and X-ray slopes 
as well as the optical-to- X-ray ratio (see Fig 16). The tem- 
perature is lower than that required for the RIXOS AGN 
but Ro.2keV is the same (see Table 4); the amount of cold 
absorption required is very small (a maximum of 3xl0 20 
cm -2 ). With regard to ultra-soft X-ray sources, we note 
that the optical-to-X-ray spectrum of the ultra-soft X-ray 
quasar E1346+266 also compares well with an unabsorbed 
bremsstrahlung spectrum (Tb rC m=7xl0 6 K; Puchnarewicz 
ct al. 1995), in agreement with our model. However for the 
lower-luminosity ultra-soft Seyfert 1 galaxy RE J1034+396, 
a bremsstrahlung plus power-law provided a poor fit to the 
data (Puchnarewicz et al. 1995), thus while this model can 



be extended to some ultra-soft X-ray AGN and UV-excess 
quasars, it may perhaps be inappropriate for lower lumi- 
nosity ultra-soft objects like RE J1034+396. However given 
the idiosyncrasies of individual objects, more stringent tests 
must be made using higher quality data before any conclu- 
sions can be drawn at this stage. 

We now discuss the implications of the model for an 
optical-to-soft X-ray big bump and for the presence and lo- 
cation of dusty, cold gas in AGN. 

5.1 The origin of the BBB 

The correlations in the data suggest that the BBB and the 
soft X-ray excess are either part of the same component (the 
'big bump') or that they are separate but subject to the 
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Figure 15. The distribution of a x versus a x for the RIXOS AGN (crosses) compared with that of the WF (diamonds) and L94 
(squares) samples. The WF and L94 data have been corrected for the differences in the definitions of « x between all three samples (see 
Section 3.1.3). 
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Figure 16. The L94 data (plotted as squares) compared with the best-fitting absorption model; (a) a x versus a x , and (b) ct op t 
versus o x . The model (plotted as a solid line) assumes a redshift of 0.2 (the mean redshift of the L94 sample) and has a temperature, 
Tbrcm=6xl0 5 K, -Ro.2fceV=10 and uses a range of JV H from to 3xl0 20 cm~ 2 . 



same external influences. In the modelling, we used a single 
optical/soft X-ray big bump and this was demonstrated to 
provide a good representation of the data, therefore through- 
out this section we have assumed a big bump as the origin 
of the BBB and the soft X-ray excess. 



5.1.1 Accretion discs 

The optically-thick AD models are a natural assumption for 
the origin of the BBB and ADs have been successfully used 
to reproduce optical-to- X-ray AGN spectra (e.g. Arnaud 
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et al. 1985; Czerny & Elvis 1987; Sun & Malkan 1989), even 
in extreme cases such as the EUV-bright Seyfert 1 galaxy 
RE J1034+396 (Puchnarewicz et al. 1995). Although some 
AD models are not appropriate for the optical/soft X-ray big 
bump in the RIXOS AGN (Section 4.1.1), we find that the 
Czerny & Elvis (1987) model of an AD with a surrounding 
hot (T — 10 8 K) corona does compare well with the most 
convex (i.e. the mean a x =1.65 spectrum indicated by the 
dashed line in Fig 9d) of the mean spectra. It is similar in 
shape to the bremsstrahlung spectrum which was success- 
fully used in the modelling. 

The AD corona is optically thin and cooled by inverse 
Compton scattering so that the high energy tail of the disc 
spectrum hardens and extends to higher energies. The ge- 
ometry of the corona is not well constrained, for example it 
may form a thin 'skin' around the disc or be large and spher- 
ical. If the optical depth of a large corona is high (r ~ 1), 
then the model's prediction that soft X-ray excesses are only 
observed from ADs observed face-on, is relaxed (Czerny & 
Elvis 1987). 

The Czerny & Elvis model shown in Fig 10 has a black 
hole mass of 10 8 Mq. With the naive assumption that the 
black hole mass scales inversely with the square root of the 
temperature (this applies in the simple thin disc approxi- 
mation, e.g. Pringle 1981), then the range of Tb rC m implied 
by the model (i.e. 10 5 — 3 x 10 6 K at z=l) suggests black 
hole masses of ~ 10 6 — 1O 9 M0 for an equivalent AD model. 
Furthermore, as the redshift of the AGN increases, the tem- 
perature of the big bump tends to decrease (see Table 4). 
In the accretion disc models, the temperature of the disc 
is determined by the mass of the black hole and is lower 
when the black hole's mass is high. Thus if the big bump 
spectrum is produced by an accretion disc, then our model 
favours more massive black holes at higher redshifts. A sim- 
ilar pattern (i.e. low-temperature ADs at high-z) would be 
expected in order for the AD model to be consistent with 
the restrictions imposed by the Eddington limit. 

5.1.2 Optically-thin gas 

The bremsstrahlung model for the big bump has many at- 
tractive features which make it a viable alternative to the 
more popular, optically-thick AD (Barvainis 1993). For ex- 
ample, it can predict observed quasar continuum slopes in 
the optical (a O pt~0.2-0.3, e.g. Neugebauer et al. 1987) with- 
out the need for an extension of the X-ray power law into the 
optical to match the observed spectra; an additional power- 
law component is often required for AD models. Also, if 
the geometry of the gas is toroidal, it can produce linearly 
polarized light similar to that which is often observed (see 
Barvainis 1993 and references therein). 

A bremsstrahlung spectrum has been assumed for the 
big bump in the modelling and we find that this success- 
fully reproduces the observed spectral indices (i.e. a ox , a opt 
and a x ). Although our model assumes that the X-ray 
power- law component continues down to ~ 10000 A, the 
bremsstrahlung spectrum dominates over the power law in 
the optical so that, in most cases, the removal of the power 
law longward of the EUV would make no significant dif- 
ference to the model fits. The optical extension to the 
power-law is required however when the total flux in the 
bremsstrahlung component is weak relative to the X-ray 



power-law (see, for example, the extreme case in Fig 14 
where there is no bremsstrahlung component). 

The scatter in the correlations observed in the RIXOS 
AGN may be satisfied by differences in the intrinsic Tb rem 
and in the relative normalization of the big bump to the un- 
derlying continuum. The best-fitting temperature of the gas 
(Tb re m~ 10 6 K) would be high enough to make the Lyman 
edge very weak and therefore difficult to detect (Barvainis 
1993), which is consistent with observations (e.g. Antonucci, 
Kinney & Ford 1989; Koratkar, Kinney & Bohlin 1992). 

One problem with an optically-thin origin of the big 
bump is its ability to reproduce the very short-timescale 
variability observed in some AGN, although Barvainis 
(1993) demonstrated that many small, free- free emitting 
cloudlets may have a source size as small as 1 lt-hr across. 

5.1.3 Reprocessing in cool clouds 

In the cool clouds model, the primary quasar spectrum is 
reprocessed by cool, dense clouds which lie close to, and 
perhaps even within, the central continuum emitting source 
(Guilbert & Rees 1988; Ferland and Rees 1988). The clouds 
absorb soft and medium X-rays and re-emit them in the 
optical/UV. This model may, for certain values of the cloud 
density and filling factor, predict observed AGN features 
such as the BBB and the soft X-ray excess, and has already 
been invoked to explain the medium-to-hard-X-ray spectra 
observed in Seyfert galaxies (Nandra & George 1994). 

Barvainis (1993) has suggested that in principle, it is 
possible for the reprocessed optical/UV continuum to be 
optically-thin, thus the reprocessing (cool) clouds them- 
selves may produce the big bump spectrum observed. The 
temperature of the free-free spectrum in this case would be 
~ 10 5 — 10 6 K, which is similar to that implied by our mod- 
elling (~ 10 5 — 3x 10 6 K), and the clouds must have relatively 
low densities (N < 10 15 cm~ 3 ; Barvainis 1993), although 
not so low that the free-free emitting source becomes very 
large. 

5.2 The absorbing medium 

The model for the RIXOS AGN optical-to-X-ray continua 
described in Section 4.1, suggests that cold gas with a 
roughly Galactic gas-to-dust ratio and moderate (A^h< 5 x 
10 21 cm~ 2 ) column density lies along the line of sight to 
the quasar nucleus. Similar levels of intrinsic absorption 
have been observed in other AGN, e.g. in the quasar 3C109, 
Allen & Fabian (1992) found evidence of cold gas with a 
column density of 5xl0 21 cm~ 2 and a Galactic dust-to-gas 
ratio, while Mathur et al. (1994) have shown that the X-ray 
and UV spectra in 3C 351 are absorbed by a highly ionized 
outflow with a column density of 1-2 xlO 22 cm~ 2 . 

Based upon the parameters of the gas suggested by the 
modelling presented in Section 4, we now explore the nature 
of a possible cold absorbing medium in the RIXOS AGN 
and whether the results of the modelling can be consistent 
with the emerging picture of quasar and Seyfert nuclei. 

5.2.1 Geometry of the gas 

A 'first-order' indication of the covering factor of this dusty 
gas may be derived from the distribution of the intrinsic Nh 
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which is implied by the modelling. 

For any given model column density (A?Hmod), the ab- 
sorbed bremsstrahlung-plus-power-law model makes a pre- 
diction of the slope of the X-ray spectrum (a xm od); a X mod 
is plotted as a function of A^mod in Fig 17a for the best- 
fitting models at redshifts from 0.2 to 2.5. Therefore, by us- 
ing the observed a x and interpolating between the points of 
the function (appropriate to the redshift of that AGN), the 
amount of absorbing material intrinsic to the AGN (A^Hint) 
may be derived. The data were thus divided into redshift 
bins (i.e. ,2=0-0.25, 0.25-0.75, 0.75-1.25, 1.25-1.75 and 1.75- 
3.5) and AT H i n t was calculated for each AGN by interpolat- 
ing between the corresponding a xmo d vs A^mod functions 
(i.e. the z=0.2, 0.5, 1, 1.5 and 2.5 curves respectively). 

The resulting distribution of A/Hint is shown in Fig 17b; 
these represent the column densities of the absorbing gas and 
dust intrinsic to the AGN and are in addition to NnGa.\- It 
shows a peak at A r H< 10 21 cm -2 and implies that approx- 
imately one-third of the RIXOS AGN lie behind columns 
> 10 21 cm" 2 . 

When considering the distribution of ATmnt, a number 
of points should be borne in mind. Firstly, the sample used 
in this paper is not fully complete, although the a x distribu- 
tion of these sources is not significantly different from that 
of the complete RIXOS survey (Mittaz et al. 1996). Sec- 
ondly, the defining criteria of the sample select against very 
soft sources which have little or no flux above 0.4 keV. How- 
ever, based on the number of AGN in the Ultra-Soft Survey 
(Thompson & Cordova 1994) which have no significant emis- 
sion above 0.5 keV, and comparing this with the AGN from 
the Einstein Extended Medium Sensitivity Survey (Gioia 
ct al. 1990; Stocke et al. 1991), we estimate that only ~10 
very soft sources have been omitted. Finally, the RIXOS 
sample will not include heavily absorbed sources which have 
insufficient 0.4-2.0 keV flux to be detected, therefore at rel- 
atively high values of A/Hint (e.g. higher than ~ 10 21 cm -2 ), 
the distribution in Fig 17b should be considered a lower 
limit. 

Taking into account these restrictions on the sample, 
we can estimate the proportion of sources which are rela- 
tively unabsorbed (i.e. with Nmnt< 10 21 cm -2 ). Assuming 
(naively) that this distribution is typical of the AGN pop- 
ulation in general, we find that 67 out of 108 sources are 
'unabsorbed' (from Fig 17), which gives a covering factor 
for the absorber of ~0.3. Allowing for a possible 10 missing 
ultra-soft X-ray sources makes no significant difference to 
the covering factor, but as we have made no allowance for 
AGN which are too strongly absorbed to be detected, this 
must be considered a very conservative lower limit. Thus 
geometries of the gas where the covering factor is relatively 
high, e.g. spherical and toroidal distributions, are favoured 
while models with a small covering factor (less than ~0.1), 
e.g. thin discs, are unlikely. 

A thick torus 

If the gas is contained within a geometrically-thick 
torus, then the range of absorbing columns implied by the 
model may be due to orientation effects if the amount of 
gas along the line of sight to the nucleus changes with the 
'viewing angle' (the angle between the line of sight and the 
axis of the torus). For example, if the density of the gas de- 
creases with increasing distance from the plane of the torus 
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Figure 17. (a) The model a x at 2=0.2, 0.5, 1.0, 1.5 and 2.5 
(solid, dotted, dashed, dot-dashed and dot-dot-dot-dashed respec- 
tively) calculated from the best-fitting absorbed bremsstrahlung 
plus power-law models at those redshifts (see Section 4.1), and 
plotted as a function of the assumed input column densities. 
(b) The distribution of inferred intrinsic column densities for 
the RIXOS AGN (iVjjint; solid line) calculated using the models 
plotted in (a) (i.e. using bins of z=0.0-0.25, 0.25-0.75, 0.75-1.25, 
1.25-1.75 and 1.75-3.5 respectively). These represent the column 
densities which would be intrinsic to the AGN themselves and are 
in addition to any Galactic absorption. 

(whether it is due to fewer clouds, less dense clouds or a 
less dense, smoothly distributed gas), then the amount of 
absorption observed would increase as the viewing angle in- 
creases. Otherwise, if the density remains constant, then 
the azimuthal profile of the torus may be such that the path 
length through the torus changes with the viewing angle. 
In this orientation-dependent model then, when the torus 
is observed relatively face-on, there would be very little ab- 
sorption and cv x is soft. As the viewing angle increases, the 
amount of cold gas along the line-of-sight increases and cv x 
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hardens. 

Alternatively, the range of Nmnt may simply indicate 
that absorbing tori in AGN have a wide range of optical 
depths. Note that if the torus is made up of individual 
clouds, then the clouds must be able to have low column den- 
sities (< 10 21 cm -2 ) in order to produce the small columns 
observed. 

Spherical 

If the gas is spherically distributed, then the observed 
range in a x reflects a range in the column density of the 
absorbing material. The distribution in Fig 17 would then 
imply that most RIXOS AGN are surrounded by gas with 
columns of ~ 10 21 cm -2 . 

5.2.2 Location of the gas 

The presence of dust places limits on the location of the 
gas, since the minimum radius at which dust grains sublime 
(r m i n ) is given by r m i n ~ 0.20L4g 5 pc, where is the bolo- 
metric luminosity in units of 10 46 erg s -1 (Laor & Drainc 
1993). For the RIXOS AGN then, the smallest radius at 
which dust can survive is the order of a parsec, which places 
the dust beyond the broad line region (BLR) . Thus the most 
likely locations of the dust are (1). the narrow line region 
(NLR), (2). the dusty molecular torus, which lies between 
the BLR and NLR, (3). the host galaxy and (4). intervening 
galaxies along the line of sight. 

NLR 

It has been suggested that dust lies within the NLR it- 
self (e.g. Laor & Draine 1993). The covering factor of the 
NLR is very small however (~0.01) and since the distribu- 
tion of iVmnt implies that a large fraction of the RIXOS 
AGN are absorbed, we conclude that the NLR is an un- 
likely source of the absorbing medium which produces the 
continuum changes. 

Molecular torus 

It has been proposed that a dusty, molecular torus lies 
between the BLR and NLR of many Seyfert galaxies (see 
e.g. Antonucci 1993 for a review). This torus is probably 
geometrically-thick [the ratio of the height (h) to the inner 
radius (r), h/r, is ~0.7, and r ~1 pc] and has a high column 
density (N H > 10 24 cm" 2 ; Krolik & Begelman 1988). It is 
made up of clouds which individually have a large column 
density (i.e. A%> 10 24 cm -2 ) and the torus has a cloud 
covering factor ~1. 

A torus of this kind is not possible as the source of 
absorbing gas in the RIXOS AGN, because the column den- 
sity of individual clouds is much higher than that implied 
by the modelling, so that observed AGN continua would 
be completely unobscured or very strongly absorbed; in 
most cases we observe moderate absorption (i.e. -/Vmnt~ 10 21 
cm -2 ). However, if the Eddington luminosity of an AGN is 
high, the torus might instead be optically thin or geomet- 
rically thin (Krolik & Begelman 1988). We have ruled out 
a geometrically-thin absorber on the basis of the A^Hint dis- 
tribution (Fig 17 and Section 5.2.1); a geometrically-thick 
torus, subject to the limitations discussed in Section 5.2.1, is 
possible however if the clouds have a low Nn (< 10 21 cm -2 ). 



Host galaxy 

The column densities implied by the modelling are typ- 
ical of the interstellar medium in the disc of our own Galaxy 
(~ 10 21 cm -2 in a direction perpendicular to the plane of 
the disc; Stark et al. 1992), however absorption by the galac- 
tic disc local to the AGN is unlikely as the covering factor 
is small. Cold gas and dust in the bulge of spiral galax- 
ies (which are the typical hosts for Seyfert 1 nuclei) and in 
ellipticals may produce the absorption implied; this would 
suggest a spherical distribution of the gas. 

Intervening galaxies 

It is also possible that the quasar spectra are absorbed 
through relatively face-on galaxies along the line of sight 
which are too faint to be seen. Spiral galaxies would be the 
most likely absorbers as they contain a higher HI mass than 
ellipticals or lenticulars (Rao & Briggs 1993). Any absorp- 
tion by external galaxies close to high-redshift quasars would 
be difficult to measure in X-rays because the strongest ef- 
fects of the absorption are shifted out of the observed range; 
the absorption we measure would be dominated by relatively 
nearby (i.e. low-z) galaxies. Nonetheless, we would expect 
to see a trend towards more heavily absorbed AGN as the 
redshift increases, yet there is no hardening of a x with z for 
the RIXOS AGN (Fig 4, this paper; Mittaz et al. 1996). 

Although it is difficult to measure this kind of absorp- 
tion in X-rays due to the redshift effect, it should be eas- 
ier in the optical because the shorter wavelength regions, 
where the dust absorption is higher, are moved into the 
observer-frame, making absorption easier to detect. In this 
case, the optical slope would soften with increasing z, yet 
again, we find no significant changes of a opt with z (for 
z > 0.25). Thus we conclude that absorption by intervening 
spiral galaxies is unlikely. 

5.3 The unified AGN model 

In the 'hidden Seyfert 1' picture of a Seyfert 2 galaxy, some 
(if not all) Seyfert 2s are actually Seyfert 1 nuclei viewed 
through the dusty molecular torus (e.g. Miller & Goodrich 
1987); this is a part of the 'unified AGN model' which essen- 
tially aims to explain the wide range of AGN properties by 
orientation effects (see e.g. Antonucci 1993 and references 
therein). We now speculate on an extension of this model 
to the quasar population in the light of the RIXOS results. 

If indeed the RIXOS AGN are absorbed as we propose, 
then what we might be seeing is a mixture of type 1 and type 
2 quasars (analogous to type 1 and type 2 Seyferts) whose 
X-ray and optical spectra are modified by a dusty torus. The 
torus has a low column density, (perhaps as a consequence 
of the relatively high luminosity of the quasars), low enough 
to modify the optical-to-X-ray continuum but still allowing 
much of the broad permitted line emission to escape. Thus 
most quasars would be classified as type Is based on their 
optical spectra, even though their nuclear emission is viewed 
through the torus (as is assumed for Seyfert 2s). This would 
explain the rarity of 'type 2' quasars (i.e. quasars with nar- 
row permitted and forbidden lines) since the BLR is not 
always completely absorbed by the torus. 

6 SUMMARY 

We have used the RIXOS AGN to probe the nature of the 
optical-to-X-ray continuum in a large sample of X-ray se- 
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lected Seyfert Is and quasars. The sample provides a wide 
range in parameter space, i.e. in luminosity, X-ray spectral 
slope and optical-to-X-ray ratio, over which to search for 
relationships between the optical and X-ray continua. We 
confirm the a ox vs q x correlation previously reported by WF 
and L94, and identify anti-correlations between a op t and a ox 
and between a op t and a x (especially at z > 1.0), suggesting 
that the optical and X-ray slopes fall and rise towards the 
EUV together. 

These relationships imply overall changes in spectral 
shape from convex (when a x is soft) to concave (when cv x 
is hard; see Fig 9d). They also suggest that the BBB and 
the soft X-ray excess are strongly dependent and they may 
be part of the same optical-to-soft X-ray 'big bump'. The 
changes in the mean optical-to-X-ray continua are not con- 
sistent with a component which shifts between the optical 
and X-ray regions, e.g. ADs or bremsstrahlung spectra with 
a range of temperatures. Although these changes might be 
produced by a big bump with varying strengths superposed 
on a relatively invariant underlying continuum, the under- 
lying continuum must be either (1) concave from 6000 A to 
2 keV (like the dot-dashed continuum in Fig 9d) yet have 
a ~1 at lower and higher frequencies, or (2) absorbed by 
cold gas and dust while the BBB and soft X-ray excess are 
completely unabsorbed: we conclude that both are unlikely. 

We have proposed instead that the changes are due to 
absorption of the entire intrinsic quasar spectrum by cold 
gas and dust and have demonstrated the viability of this 
idea from 2=0 to ~3 by creating models of the absorbed 
continua. The model assumes that the intrinsic spectrum is 
the sum of an a=l power law and a bremsstrahlung spec- 
trum. The temperature of the bremsstrahlung component 
favoured by the modelling is ~ 10 K and the scatter in 
the data implies a range in Tb rC m of ~ 10 5 K to 3x 10 6 K. A 
wide range in the strength of the big bump is also permitted 
by the model, e.g. at 2=1 and for T brom =1.25xl0 6 K, the 
ratio of flux in the bremsstrahlung spectrum to that in the 
power-law at 0.2 keV can vary from to ~1000. Thus the 
model does not rule out the possibility of significant differ- 
ences in the strength and temperature of the big bumps from 
source to source. It does suggest however that the correla- 
tions observed are not due to changing big bump parameters 
(i.e. neither the strength nor temperature) but that they are 
consistent with the effects of absorption by dusty cold gas. 

The big bump may be produced by an AD with a sur- 
rounding hot corona, an optically-thin gas or perhaps re- 
processing in cool clouds. The absorber is cold gas with an 
approximately Galactic dust abundance and its column den- 
sity varies from to ~ 4 x 10 21 cm~ 2 . The distribution of the 
intrinsic Nh in the RIXOS AGN implied by the modelling 
suggests a high covering factor, e.g. a toroidal or spherical 
geometry of the gas. The presence of dust places it beyond 
the BLR and, if it is toroidal, then the different types of 
observed continua may be due to an orientation effect with 
the column density increasing with the viewing angle. 

If the obscuring torus which lies between the BLR and 
the NLR is responsible for the absorption, then the column 
density of the clouds which make up the torus must be signif- 
icantly lower than that observed in Seyfert 2s. We speculate 
that the rarity of 'type 2' quasars is due to obscuring tori 
with low column densities which do not always fully absorb 
the BLR; this may be a physical consequence of the high 



luminosity in quasars compared to Seyferts. An alternative 
source of the absorbing gas may be the bulge of a spiral host 
galaxy or an elliptical host galaxy. 
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